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Frequency  Dependent  Shear-wave  Q  Models 
for  the  Crust  of  China  and  Surrounding  Regions 

Alemayehu  L.  Jemberie  and  Brian  J.  Mitchell 
Department  of  Earth  and  Atmospheric  Sciences 
Saint  Louis  University 

Abstract 

Starting  with  fundamental-mode  Rayleigh-wave  attenuation  coefficient  values  (yr) 
predicted  by  previously  determined  frequency-independent  models  of  shear-wave  Q  ( Q ),) 
we  have  obtained  frequency-dependent  Q  models  that  explain  both  the  measured  values 
of  yr  as  well  as  of  Lg  coda  Q  and  its  frequency  dependence  at  1  Hz  ( Q0  and  ? 7 
respectively)  for  China  and  some  adjacent  regions.  The  process  combines  trial-and-error 
selection  of  a  model  for  the  depth  distribution  of  the  frequency  dependence  parameter  (£) 
for  Q  with  a  formal  inversion  for  the  depth  distribution  of  Q ^  at  1  Hz.  Fifteen  of  the 
derived  models  have  depth  distributions  of  £  that  are  constant,  or  nearly  constant, 
between  the  surface  and  a  depth  of  30  km.  ^distributions  that  vary  with  depth  are 
necessary  to  explain  the  remaining  seven  models.  £  values  for  the  depth-independent 
models  vary  between  0.4  and  0.7  everywhere  except  in  the  western  portion  of  the  Tibetan 
Plateau  where  they  range  between  0.1  and  0.3  for  three  paths.  These  low  £  values  lie  in  a 
region  where  QLg  and  crustal  Qfl  are  very  low  and  suggest  that  they  should  also  be  low  for 
high-frequency  propagation.  The  models  in  which  £  varies  with  depth  all  show  a  decrease 
in  that  value  ranging  between  0.6  and  0.8  in  the  upper  15  km  of  the  crust  and  (with  one 
exception  where  £  =  0.0)  between  0.3  and  0.55  in  the  depth  range  15-30  km.  The 
distribution  of  highest  £  values  (0.6  -  0.8)  in  the  upper  crust  indicates  that  high-frequency 


2 


waves  will  propagate  most  efficiently,  relative  to  low-frequency  waves,  in  a  band  that 
includes,  and  strikes  north-northeastward  from,  the  path  between  event  212/97  and  KMI 
to  the  path  between  event  1 80/95  and  station  HIA  in  the  north. 
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Introduction 

Almost  all  mechanisms  for  intrinsic  Q  in  the  Earth  predict  that  it  should  vary  with 
frequency  (e.g.  Jackson  and  Anderson,  1970).  The  variation  of  Q  with  frequency  is 
usually  considered  to  vary  as  Q  ~  ooa  where  the  exponent  a  is  termed  the  frequency 
dependence  parameter  and  may  be  different  for  different  attenuating  waves.  In  the 
present  study  we  represent  that  parameter  for  shear-wave  Q  (Qfl)  by  t  and  for  Lg  coda  Q 
(QcLg)  by  r|.  Several  studies  have  addressed  the  question  of  frequency  dependence  of  Q  in 

the  mantle  and  have  found  it  to  be  necessary  for  simultaneously  explaining  the 
attenuation  rate  of  both  low-  and  high-frequency  waves.  Studies  have  invoked  frequency 
dependence  to  reconcile  Q  values  observed  for  free  oscillations  of  the  Earth  and  1-Hz 
body  waves  (Jeffreys,  1967;  Liu  et  al.,  1976),  to  explain  shear-wave  Q  observed  for  ScS 
waves  at  low  and  high  frequencies  (Sipkin  and  Jordan,  1979),  and  to  explain  variations  in 
Q  with  frequency  for  teleseismic  P-  and  S-waves  in  the  range  0.02-4.0  Hz  (Der  et  al., 
1982).  Frequency  dependence  of  Q  for  upper  mantle  material  rock  has  also  been 
observed  in  laboratory  experiments  (e.g.  Gueguen  et  al.,  1989;  Jackson  et  al.,  2003). 

Fewer  studies  have  addressed  the  question  of  frequency  dependence  for  Q  at  crustal 
depths.  Combined  inversions  of  surface-wave  attenuation  at  periods  between  about  5  and 
50  s  and  of  Lg  coda  Q  at  1  Hz,  have  shown  that  Q  in  the  continental  crust  varies  with 
frequency  at  least  at  frequencies  of  about  1  Hz  and  above  (Mitchell,  1980;  Mitchell  and 
Xie,  1994).  Determinations  of  Lg  Q  (Qig)  (e.g.  Benz  et  al.,  1997)  or  Lg  coda  Q  (QcLg) 

(e.g.  Xie  and  Mitchell,  1990a,  1990b)  typically  suggest  the  need  for  frequency-dependent 
values.  It,  however,  has  also  been  found  possible  to  explain  frequency-dependent  QLg  or 
QcLg  with  a  frequent-independent  Q ^  model  that  contains  a  rapid  increase  in  value  at  mid- 
crustal  depths  (Mitchell,  1991). 

Theoretical  work  has  shown  that  Q  frequency  dependence  may  be  associated  with 
intrinsic  energy  loss  due  to  temperature  (e.g.  Gueguen,  1 989)  or  to  movements  of  fluid  in 
permeable  rock  (e.g.  O’Connell  and  Budiansky,  1977).  Recent  studies  (Hammond  and 
Humphreys,  2000;  Paul  et  al.,  2003)  have  concluded  that  the  latter  mechanism  is  not 
plausible  at  mantle  depths  where  pressures  are  high,  but  research  on  this  question,  to  our 
knowledge,  has  not  yet  been  conducted  for  crustal  rock  at  lower  pressures.  Some 


4 


calculations,  using  realistic  values  for  temperature  and  the  frequency  dependence 
parameter  for  Q,  as  well  as  assumed  values  for  Q  activation  energy,  suggest  that 
temperature  differences  cannot  explain  regional  differences  between  crustal  Q  values  for 
the  eastern  and  western  United  States  (Mitchell,  1995).  If  this  result  is  correct  it  suggests 
that  regional  Q  variations  in  the  upper  crust  are  likely  to  be  caused  by  regional  variations 
in  fluid  content  in  crustal  rock.  Frequency  dependence  of  Q,  at  least  at  high  frequencies, 
may  also  be  attributed  to  scattering  from  heterogeneities  in  material  traversed  by  seismic 
waves  (e.g.  Dainty  et  al.,  1987). 

That  Q  varies  with  frequency,  at  least  in  some  portions  of  he  crust  and  upper  mantle 
and  over  the  frequency  range  that  that  relevant  to  seismic  wave  propagation  over  near  and 
regional  distances,  appears  to  be  well  established.  Consequently  this  frequency 
dependence  may  be  of  great  practical  consequence  for  magnitude  determinations  at  near 
and  regional  distances;  and  in  nuclear  test  ban  treaty  monitoring,  for  determining 
detection  thresholds  and  for  the  implementation  of  methods  for  discriminating  between 
earthquakes  and  explosions. 

In  this  study  we  determine  frequency  dependence  values  (£)  for  crustal  Qjt  that 
simultaneously  explain  the  attenuation  observed  for  fundamental-mode  Rayleigh  waves 
at  intermediate  periods  (about  5-50  s)  and  values  of  Unobserved  in  the  southeastern 

portion  of  Eurasia  that  includes  China  and  some  adjacent  regions.  This  will  allow  us  to 
extend  the  applicability  of  previously  determined  frequency-independent  Qfl  models 
determined  from  surface  waves  to  higher  frequencies  that  characterize  waves  recorded  at 
near  and  regional  distances. 

Our  approach  differs  from  that  used  in  previous  surface-wave  studies  of  Qfl 
frequency  dependence  in  that  we  were  able,  for  each  path,  to  use  fundamental-mode 
Rayleigh-wave  data  and  higher- frequency  Lg  coda  data  recorded  by  the  same  instrument. 
Earlier  studies,  before  the  widespread  availability  of  broadband  data  (e.g.  Mitchell,  1980) 
used  different  instruments  for  low-  and  high-frequency  data.  The  top  trace  in  Figure  1 
shows  ground  motion  recorded  by  the  broadband  station  ULN  in  Mongolia  for  an 
earthquake  about  530  km  distant  near  Lake  Baikal.  A  seismogram,  if  recorded  by  a  long- 
period  instrument  of  the  World-Wide  Standard  Seismograph  Stations  (WWSSN),  would 
be  similar  to  the  middle  trace  of  Figure  1  which  is  obtained  by  low-pass  filtering  the 
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upper  trace  at  periods  above  4  s.  If  recorded  by  a  short-period  instrument  of  the  WWSSN 
the  seismogram  would  look  like  the  bottom  trace  of  Figure  1  which  was  obtained  by 
band-pass  fdtering  the  upper  trace  around  1  Hz. 

Determination  of  Q  structure  using  fundamental-mode  Rayleigh  waves 

Using  a  single-station  multimode  method,  Jemberie  and  Mitchell  (2003)  determined 
three-layer  frequency-independent  Qft  models  for  the  southeastern  portion  of  Asia  that 
includes  China  and  some  adjacent  regions.  It  used  surface-waves  that  were  well-recorded 
at  surface- wave  frequencies  (about  5-50  s),  an  example  of  which  appears  in  the  middle 
trace  of  Figure  1.  They  mapped  Qf,  variations  for  two  depth  ranges,  0-10  and  10-30  km 
and  also  estimated  in  the  depth  range  30-60  km  but  with  less  reliability.  They  found 
large  regional  variations  of  QfX  in  all  of  the  layers.  Qfl  in  layer  1  (the  upper  10  km)  varies 
between  about  250  in  southeastern  China  to  about  40  in  thoughout  most  of  the  western 
Tibetan  Plateau.  In  layer  2  (10  -  30  km  depth)  Qu  varies  between  about  140  in  central  and 
eastern  China  and  about  60  in  western  Tibet  and  the  Burma-Thailand  region.  The  low  Q 
values  in  western  Tibet  are  among  the  lowest  found  in  any  continental  region  and 
correlate  spatially  with  the  lowest  values  of  Qig  ever  to  be  reported  (Xie,  2002). 

In  order  to  obtain  frequency-dependent  Qfl  models  our  process  requires  that  we 
invert  curves  of  fundamental-mode  Rayleigh-wave  attenuation  coefficient  values  (y) 
versus  period.  We  did  not  measure  these  directly,  but  instead  computed  y  values  for  the 
Qlt  models  of  Jemberie  and  Mitchell  (2003)  along  with  appropriate  velocity  models.  This 
computation  is  the  first  step  in  the  inversion  process  described  in  a  later  section. 

Lg  Coda  Q 

Lg  is  usually  the  most  prominent  high-frequency  phase  (-0.5  Hz  and  higher) 
recorded  on  seismograms  at  regional  distances  in  continents  and  can  be  observed  at 
distances  as  great  as  3000  km.  Because  it  is  so  well  recorded  it  has  been  used  to 
determine  magnitudes  of  small  events  at  regional  distances  (Nuttli,  1973)  and  to  study 
regional  variations  of  attenuative  properties  of  the  crust.  In  stable  regions  it  is 
characterized  by  group  velocities  of  about  3.5  km/s  and  can  be  as  low  as  3.2  or  3.3  km/s 
in  tectonically  active  regions. 


Lg  can  be  represented  as  a  superposition  of  many  higher-mode  Rayleigh  waves  at 
high  frequencies.  Computed  Lg  has  a  relatively  sharp  onset  followed  by  oscillations,  or 
coda,  that  can  continue  for  several  minutes  if  the  crustal  model  used  for  the  computations 
contains  low-velocity  layers  in  the  uppermost  crust  (see  the  bottom  trace  in  Figure  1). 
The  coda  of  observed  Lg  waves  usually  continues  for  a  longer  time  than  can  be 
theoretically  predicted  by  theoretical  seismogram  computations  using  plane-layer  models. 
These  additional  oscillations  are  generally  attributed  to  scattering  of  wave  energy  from 
crustal  heterogeneities.  Paradoxically,  even  though  the  coda  of  seismic  waves  is  largely 
due  to  scattering,  theoretical  and  computational  work  indicates  that  measurements  of  Q 
using  coda  yield  values  of  intrinsic  Q  rather  than  scattering  Q  (e.g.  Wennerberg,  1993). 
For  this  reason  we  combine  surface-wave  attenuation  at  intermediate  periods  (with 
wavelengths  likely  to  be  much  longer  than  dimensions  of  heterogeneities)  with  Lg  coda, 
at  frequencies  (near  1  FIz)  that  are  likely  to  be  affected  by  scattering,  to  invert  for 
frequency-dependent  models  of  Qfl  for  the  crust. 

Our  method  for  determining  crustal  Qfl  models  is  described  in  the  following  section. 
Those  models  should  explain  measured  values  for  fundamental-mode  Rayleigh-wave 
attenuation  as  well  as  average  values  of  Q0  and  i]  for  QcLg  along  the  paths  of  propagation. 

The  utilization  of  Rayleigh-wave  attenuation  was  discussed  in  the  previous  section. 
Tomographic  maps  of  Q0  and  r]  at  1  Hz  are  available  for  almost  all  of  Eurasia,  including 
our  region  of  study  (Mitchell  et  al.,  197).  The  Qc  values  vary  between  about  150  and 
1000  with  the  lower  values  occurring  in  a  broad  band,  the  Tethysides  belt,  that  stretches 
across  southern  Asia,  including  China  and  some  adjacent  regions  pertinent  to  this  study. 

Method 

As  a  first  step  toward  determining  the  frequency  dependence  of  Q ^  and  how  it  may 
vary  with  depth  in  China  .we  have  developed  a  method  by  which  Q j,  frequency 
dependence  can  be  obtained  for  our  crustal  models  using  a  combined  forward 
modeling/formal  inversion  procedure.  It  is  a  variation  of  a  previously  developed 
procedure  (Mitchell  and  Xie,  1994)  and  obtains  a  depth  distributon  of  £  that  reconciles  Qu 
values  at  1  Hz  with  those  obtained  at  surface-wave  periods  (about  5-50  s).  Because 
simultaneous  inversions  for  Qh  and  £  include  more  parameters  than  do  inversions  for  Q ^ 
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alone,  an  additional  level  of  non-uniqueness  is  introduced  in  the  models.  The  process 
proceeds  as  follows: 

1 .  Compute  attenuation  coefficients  for  the  fundamental  Rayleigh  mode,  as  a  function  of 
period,  for  the  pertinent  frequency-independent  Qft  models  Jemberie  and  Mitchell 
(2003)  and  the  corresponding  velocity  models.  For  our  region  of  study  we  use  the 
shear-velocity  models  from  Jemberie  (2002)  except  for  Tibet  where  we  used  a  model 
from  Chun  and  Yoshi  (1977). 

2.  Assume  a  depth  distribution  for  the  frequency  dependence,  £,  of  Q)t  and  use  the 
attenuation  coefficients  obtained  in  step  1  to  obtain  a  Qfl  model. 

3.  Calculate  the  fundamental-mode  Rayleigh-wave  attenuation  coefficient  values  which 
are  predicted  by  the  model  obtained  in  step  2  and  compare  them  with  the  values  of 
step  1.  An  example  appears  in  Figure  2.  If  they  agree  within  reasonable  uncertainties 
go  to  step  4;  if  not,  try  a  new  distribution  of  £. 

4.  From  the  results  of  step  3,  compute  Qfl  at  frequencies  of  interest  (1  Hz  for  the  present 
study).  Figure  3  shows  an  example  of  1Hz  Qtl  models  computed  using  five  different  £ 
distributions,  assumed  to  be  constant  with  depth,  for  a  path  between  station  BRVK  and 
an  event  that  occurred  on  day  325  of  the  year  1988  (325/98  in  Figures  5  and  6). 

5.  Compute  synthetic  Lg  seismograms  (Figure  4)  at  several  distances  from  a  seismic 
source  using  the  appropriate  velocity  and  Q  models  from  step  2. 

6.  Apply  the  stacked  spectral  ratio  (SSR)  method  of  Xie  and  Nuttli  (1988)  to  the  set  of 
synthetic  seismograms  in  step  4  to  obtain  values  of  Qa  and  77  predicted  by  the  derived 
Qt,  model  and  compare  them'  to  the  measured  values  from  maps  of  (Mitchell  et  al., 
1997)  in  Figures  5  and  6. 

7.  If  necessary,  change  the  depth  distribution  of  £  in  step  2  and  repeat  the  procedure. 

Table  1  shows  derived  Qa  and  77  values  for  five  different  £  values  for  a  path  between 
event  325/98  and  station  BRVK.  Values  of  both  Qa  and  77  increase  as  the  value  of  £ 
increases.  Comparison  of  Q0  values  in  Table  1  with  Figure  5  shows  that  the  best  value  of 
£  that  can  predict  Qa  lies  between  0.4  and  0.6.  We  found  (Figure  7d)  that  a  £  value  of 
0.55  best  predicts  the  value  of  Qa  observed  for  that  path  (647  ±51)  and  also  predicts  the 
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r)  value  of  0.37  ±  0.03  This  value  is  lower  than  the  average  £  value  (about  0.5  on  average) 
for  the  path  in  Figure  6.  In  cases,  such  as  this,  where  Q0  and  77  cannot  both  be  fit 
precisely,  we  assume  that  Q0  is  more  likely  to  be  measured  more  accurately  than  77.  77 
values  may  be  inaccurate  because  they  are  obtained  by  differencing  values  of  Q  at  two 
different  frequencies. 


Results 

The  frequency-dependent  models  obtained  using  the  above  procedure  appear  in 
Figures  7a-d.  The  solid  lines  on  the  left-hand  side  of  each  pair  of  boxes  are  frequency- 
independent  Qt,  models  obtained  using  the  single-station  multimode  method  (Jemberie 
and  Mitchell,  2003)  and  the  dashed  lines  show  models  at  1  Hz  using  our  procedure  for 
obtaining  Qfi  frequency  dependence.  The  boxes  on  the  right-hand  side  of  each  pair  show 
the  assumed  distributions  of  £  with  depth.  For  15  of  the  22  models,  a  constant  value  of  £, 
or  a  distribution  that  varies  with  depth  by  no  more  than  0.05,  produced  acceptable  results. 
Values  for  these  models  range  between  0.4  and  0.7  everywhere  in  the  region  of  study 
except  the  western  portion  of  the  Tibetan  Plateau  where  values  lie  in  the  range  0.1  -  0.3. 
One  path,  284/94  to  LSA  in  the  eastern  portion  of  the  Tibetan  Plateau  has  a  R  value  0.45 
which  is  intermediate  between  values  for  paths  in  the  western  part  of  the  plateau  and 
paths  outside,  and  to  the  east  of,  the  plateau. 

The  remaining  seven  models  require  that  £  values  at  depths  greater  than  1 5  km  be 
smaller  those  in  the  upper  15  km  of  the  cmst.  For  the  upper  15-km  thick  layer  £  ranges 
between  0.6  and  0.8,  and  for  greater  depths  ranges  (with  one  exception)  between  0.3  and 
0.55.  The  exceptional  value  is  0.0.  Values  of  Q0  and  77  for  Lg  attenuation  at  1  Hz 
predicted  by  the  models  appear  with  each  pair  of  panels.  Note  that,  as  shown  by  the 
model  in  the  upper  left  of  Figure  7a,  a  model  with  a  nearly  frequency-independent  can 
still  produce  a  frequency-dependent  Qig  value  if  Qft  increases  rapidly  with  depth  at  mid- 
crustal  depths.  This  phenomenon  was  shown  to  occur  in  the  crust  of  the  Basin  and  Range 
province  in  the  western  United  States  (Mitchell,  1 99 1 ): 

The  low  values  for  £  in  the  western  Tibetan  Plateau  coincide  with  very  low  values  of 
found  in  the  crust  of  this  region  using  multimode  surface  waves  (Jemberie  and 
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Mitchell,  2003),  as  well  as  with  very  low  QLg  values  at  1-Hz  frequencies  (Xie,  2003). 
These  results  indicate  that  surface-waves  at  periods  of  -5-50  s  as  well  as  Lg  or  Lg  coda 
at  1  Hz,  and  probably  higher  frequencies,  propagate  less  efficiently  in  western  Tibet  than 
anywhere  else  in  our  region  of  study.  The  rates  of  attenuation  in  this  region  for  both 
surface-wave  periods  (-5-50  s)  and  Lg  frequencies  (-about  1  Hz)  are  the  highest  so  far 
measured  for  any  continental  region.  An  early  study  of  Q  in  the  central  Tibetan  Plateau 
(Chang  and  Yao,  1979),  using  spectral  rations  of  S  and  ScS  waves,  found  a  Q  value  of 
about  25  that  must  pertain  predominantly  to  upper  mantle  depths. 

An  interpretation  of  our  results  in  terms  of  a  continuous  relaxation  model  for  the 
crust  would  indicate  that  the  relaxation  spectrum  ((Y1)  is  shifted  to  higher  frequencies 

compared  to  other  parts  of  our  region  of  study,  an  effect  that  could  be  produced  by  higher 
temperatures  or  elevated  levels  of  tectonic  stress.  Sublithospheric  mantle  heat  flow 
estimated  by  Artemieva  and  Mooney  (2001)  is  elevated  throughout  the  western  and 
central  parts  of  the  Tibetan  Plateau  relative  to  surrounding  regions.  It  is  probably 
pertinent  that  electrical  resistivity  values  at  mid-  and  lower-crustal  depths  are  very  low 
there  (Wei  et'al.,  2001),  suggesting  the  presence  of  melt  or  fluids  that  may  been  released 
by  hydrothermal  reactions  initiated  by  the  high  temperatures. 

Conclusions 

We  have  obtained  frequency-dependent  models  of  Q fl  that  explain  the  variation  of 
fundamental-mode  Rayleigh-wave  attenuation  coefficients  with  period  predicted  by  the 
Qp  models  of  Jemberie  and  Mitchell  (2003)  as  well  as  previously  mapped  values  of  Q0 
and  ?7  in  China  and  some  adjacent  regions.  £  values  for  depth-independent  models  vary 
between  0.4  and  0.7  through  most  of  the  region  of  study,  but  range  between  0. 1  and  0.3  in 
the  western  portion  of  the  Tibetan  Plateau.  The  latter  (low)  £  values  coincide  with  regions 
where  Qig  and  crustal  Q ^  were  previously  reported  to  be  very  low  and  indicate  that  high- 
frequency  propagation  should  also  be  low  compared  to  other  regions.  The  models  in 
which  £  varies  with  depth  all  show  a  decrease  in  that  value  ranging  between  0.6  and  0.8 
in  the  upper  15  km  of  the  crust  and  between  0.3  and  0.55  (with  one  exceptionally  value) 
in  the  depth  range  15-30  km.  The  distribution  of  highest  £  values  (0.6  -  0.8)  in  the  upper 
crust  indicates  that  high-frequency  waves  will  propagate  most  efficiently,  relative  to  low- 
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frequency  waves,  in  a  band  that  includes,  and  strikes  north-northeastward  from,  the  path 
between  event  212/97  and  KM1  to  the  path  between  event  180/95  and  station  HI  A  in  the 
north. 

Our  results  indicate  that  wave  propagation  at  surface-wave  periods  (-5-50  s),  at  1 
Hz,  and  probably  higher  frequencies,  are  less  efficient  in  western  Tibet  than  anywhere 
else  in  our  region  of  study  and  that  the  region  is  characterized  by  the  least  efficient 
propagation  of  seismic  waves,  at  all  seismic  frequencies,  yet  known  for  any  continental 
region.  The  low  t,  values  can  be  explained  if  wave  propagation  in  the  caist  of  this  region 
can  be  characterized  by  a  continuous  relaxation  model  that  has  been  displaced  to  higher 
frequencies  by  temperatures  and/or  tectonic  stress  levels  that  are  higher  than  they  are  in 
surrounding  regions. 
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Table  Caption 

Table  1.  Lg  Q  at  1  Hz  ( Qo )  and  frequency  dependence  ( rj )  values  derived  from  the 
stacked  spectral  ratios  (Xie  and  Nuttli,  1988)  of  the  synthetic  seismograms  obtained 
using  the  1-Hz  shear-wave  Q  ((?,,)  structure  and  its  frequency  dependence  (t)  for  the 
path  between  event  325/98  and  station  BRVK.  The  numbers  in  column  for  F  give  the 
frequency  ranges  utilized  in  the  determinations. 

Figure  Captions 

Figure  1.  Vertical-component  seismogram  recorded  by  the  broadband  instrument  at 
station  ULN  for  an  earthquake  that  occurred  on  30  June  1995  at  an  epicentral 
distance  of  about  530  km.  The  top  trace  shows  the  original  broadband  ground 
motion,  the  middle  trace  is  the  original  record  after  it  has  been  low-pass  filtered  for 
frequencies  below  0.25  Hz  (periods  above  4  s),  and  the  bottom  trace  is  the  original 
record  after  it  has  been  band-pass  filtered  around  1  Hz. 

Figure  2.  Comparison  of  Rayleigh-wave  attenuation  coefficients  predicted  by  the  three- 
layer  frequency-independent  model  of  with  a  frequency-dependent  model  for 
spectra  obtained  at  station  ULN  for  event  180/95.  The  necessary  t  values  are  0.65 
for  depths  between  0  and  15  km  and  0.60  for  depths  greater  than  15  km. 

Figure  3.  1-Hz  models  computed  using  five  different  values  of  Qf,  frequency 
dependence  (£). 

Figure  4.  Synthetic  Lg  seismograms  generated  using  the  Qfl  structure  obtained  at  1  Hz 
and  the  velocity  model  for  a  subregion  of  China. 

Figure  5.  Topographic  map  of  Lg  coda  Q  at  1  Hz  (Q0)  for  most  Eurasia  (modified  from 
Mitchell  et  al.,  1997)  showing  paths  between  events  and  stations  where  frequency- 
independent  models  of  Qfl  structure  were  determined  by  Jemberie  and  Mitchell 
(2003). 

Figure  6.  Tomographic  map  of  the  frequency  dependence  (q)  of  Lg  coda  Q  at  1  Hz  for 
most  of  Eurasia  (modified  from  Mitchell  et  ah,  1997)  showing  paths  between  events 
and  stations. 

Figure  7a.  Frequency-dependent  models  of  Q  structure  for  Tibet  and  the  southernmost 
portion  of  the  study  area.  The  solid  line  in  left-hand  panels  for  each  path  presents 
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frequency-independent  Q  models  that  explain  measured  Rayleigh-wave  attenuation 
at  periods  between  about  5  and  50  s.  The  model  delineated  by  the  dashed  line  shows 
a  Q  model  at  1  Hz  that  explains  both  fundamental-mode  Rayleigh-wave  attenuation 
and  measured  Lg  coda  Q  values  at  1  Hz.  The  caption  identifies  the  recording  station 
and  date  of  the  event  and  gives  values  for  Lg  Q  and  its  frequency  dependence  at  1 
Hz  that  are  predicted  by  the  1-Hz  Q  model.  The  right-hand  panel  shows  the  depth 
distribution  of  Q  frequency-dependence  that  explains  the  two  models  on  the  left. 
The  caption  in  this  panel  gives  the  value  of  £,  for  the  upper  15  km  of  the  model  (first 
number)  and  for  the  depth  range  between  15  and  60  km  (second  number). 

Figure  7b.  Frequency-dependent  models  of  Qft  structure  for  central  and  western  China, 
Mongolia  and  peripheral  regions.  See  the  caption  for  Figure  7a  for  an  explanation  of 
the  models  and  symbols. 

Figure  7c.  Frequency-dependent  models  of  Qp  structure  for  paths  in  northern  China, 
Mongolia  and  portions  of  southern  Siberia.  See  the  caption  for  Figure  7a  for  an 
explanation  of  the  models  and  symbols. 

Figure  7d.  Frequency-dependent  models  of  Qp  structure  for  the  western  portion  of  the 
study  area.  See  the  caption  for  Figure  7a  for  an  explanation  of  the  models  and 
symbols. 
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Computational  Study  of  High-frequency  P-wave  Synthetics 
and  their  Attenuation  in  Continental  Crust 

F.  Leyton,  R.  Chu,  A.  Fatehi  and  B.J.  Mitchell 
Department  of  Earth  and  Atmospheric  Sciences 
Saint  Louis  University 

Abstract 

We  have  completed  preliminary  studies  of  the  effect  that  simple  velocity  and  Q 
models  have  on  the  attenuation  of  high-frequency  P-waves  for  distances  between  10  and 
200  km  from  a  source  in  the  upper  crust.  We  computed  complete  waveforms  that  might 
propagate  in  simple  high-  and  low-Q  crustal  models,  one  corresponding  to  a  typical  stable 
region  and  the  other  to  a  low-Q  model  that  approximates  a  portion  of  a  region  of  plate 
collision  in  China.  Both  models  consist  of  a  10-km  thick  upper  crustal  layer  over  a  half¬ 
space.  Amplitudes  were  determined  as  a  function  of  distance  by  computing  record 
sections  of  complete  waveforms  using  a  wave-number  integration  code.  The  computed 
seismograms  show  that  the  rate  of  attenuation  with  distance  is  much  more  rapid  for  the 
low-Q  model  than  for  the  high-Q  model  at  the  higher  frequencies  but  only  slightly  more 
rapid  at  1-Hz.  Record  sections  show  that  the  maximum  amplitudes  occur  at  epicentra] 
distances  somewhat  smaller  that  the  critical  distance  for  the  model  we  used.  Maximum 
amplitudes  for  waveforms  that  have  been  band-pass  fdtered  at  1,  3  and  5  Hz  shift  to 
lower  frequencies  with  decreasing  distance.  Our  computations  show  that  the  value 
selected  for  the  frequency-dependence  parameter  (£)  will  have  a  very  large  effect  on  the 
computed  amplitudes.  For  vertical-component  ground  motion  our  selected  t,  value  (0.5) 
yields  maximum  amplitudes  at  5  Hz  that  are  about  4  times  as  large  as  those  at  1  Hz  and 
1 .5  times  as  large  as  those  at  3  Hz.  Computations  of  horizontal-component  ground 
motion  yield  similar  results. 
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Introduction 

Studies  of  seismic-wave  attenuation  have,  over  that  past  20-30  years,  greatly 
enhanced  our  understanding  of  the  anelastic  properties  of  the  crust  and  upper  mantle  at 
surface-wave  frequencies  (-0.01-0.2  Hz)  as  well  as  for  some  waves  at  higher  frequencies, 
such  as  Lg  and  its  coda  near  1  Hz.  Much  less  is  known  about  the  attenuation  of  seismic 
waves  frequencies  higher  than  1  Hz.  These  waves  have,  however,  taken  on  greater 
importance  as  test-ban  monitoring  programs  attempt  to  use  shorter  paths  than  ever  before 
in  order  to  detect  and  locate  small  events  and  to  distinguish  them  from  earthquakes. 

We  need  to  study  the  attenuation  of  these  phases  if  we  are  to  understand  all  facets  of 
Q  and  its  frequency  dependence  over  broad  frequency  ranges.  The  purpose  of  this  section 
of  our  report  is  to  describe  preliminary  results  obtained  from  the  computation  of  P-wave 
synthetics  for  simple  velocity  and  Q  models  that  might  be  characterized  by  reasonable 
values  for  velocity  and  Q  at  relatively  high  frequencies  (1-10  Hz). 

We  computed  waveforms  for  direct  P-waves  using  a  wave-number  integration  code 
developed  at  Saint  Louis  University  and  MIT  by  B.  Mandal  (Mandal  and  Mitchell,  1986; 
Mandal  and  Toksoz,  1990).  While  the  early  codes  synthesized  complete  waveforms  for 
the  specialized  case  of  transverse  isotropy  with  a  vertical  axis  of  symmetry,  the  more 
recent  codes  perform  those  computations  for  waves  traversing  generalized  (up  to  21 
elastic  constants)  anisotropic  media.  While  we  have  performed  several  types  of 
computation,  some  with  the  intention  of  investigating  the  possibility  of  anisotropic  Q  in 
the  continental  crust,  the  present  application  restricts  computations  to  wave  propagation 
in  isotropic  media. 

We  consider  our  computations  to  be  preliminary  but  they  have  yielded  results 
suggesting  that  knowledge  of  high-frequency  wave  attenuation  will  be  an  important  area 
for  future  research  directed  toward  monitoring  a  comprehensive  nuclear-test-ban  treaty  or 
CTBT. 

High-frequency  P-wave  Attenuation  for  High-Q  and  Low-Q  Crustal  Models 

In  order  to  try  to  understand  the  attenuation  of  high-frequency  (above  1  Hz)  P-waves 
we  restricted  our  computations  in  these  preliminary  computations  to  direct  waves  at 
relatively  short  distances.  We  constructed  record  sections  for  direct  P-waves  generated  by 


a  5-km  deep  source  and  that  have  traveled  200  km  through  simple  two-layer  models 
(Table  1).  The  upper  layer  is  10  km  thick  and  has  a  P-velocity  typical  of  granitic  rock, 
with  corresponding  values  for  S-velocity  and  density. 

For  our  computations  we  took  two  end-member  cases,  a  high-Q  model  that  might  be 
representative  of  an  old  stable  continental  region  and  a  low-Q  model  that  might 
characterize  crust  in  a  zone  of  continent-continent  collision  such  as  that  occurring  in  the 
Tethysides  belt  of  southern  Eurasia.  We  have  based  the  low-Q  model  on  values  that 
might  be  expected  for  a  portion  of  southern  Asia  in  which  Jemberie  and  Mitchell  (2003) 
studied  shear-wave  Q  variations  using  Rayleigh  waves.  For  our  computations  we  assume 
that  P-wave  Q  is  twice  as  high  as  S-wave  Q  for  those  models.  These  considerations  give  a 
high-Q  P-wave  model  for  which  Q  is  1000  in  the  upper  10  km  of  the  model  and  600  in 
the  lower  layer  and  a  low-Q  model  with  a  P-wave  Q  value  of  100  in  the  upper  layer  and 
140  in  the  lower  layer. 

For  all  computations  we  assume  a  value  for  the  frequency  dependence  parameter  of 
Q  (see  part  I  of  this  report)  to  be  0.5.  This  value  is  a  rough  average  of  the  upper  crustal 
frequency-dependence  parameters  found  by  Jemberie  and  Mitchell  (2003b)  for  the  crust 
of  southeastern  Asia  that  includes  China  and  surrounding  regions.  Our  computational 
results  are  likely  to  depend  strongly  on  the  value  of  this  parameter,  so  it  will  be  important 
to  determine  its  value  and  regional  variation  for  frequency  ranges  that  are  important  in 
test  ban  monitoring. 

After  computing  the  P  waveform  at  each  distance,  we  filtered  the  trace  with  narrow 
band-pass  filters  centered  at  l  Hz,  3  Hz  and  5  Hz.  Record  sections  for  these  three  cases 
appear  in  Figures  1,  2  and  3,  each  showing  seismograms  for  both  the  high-Q  and  low-Q 
models.  The  record  sections  cover  the  distance  range  10-200  km  with  a  station  spacing  of 
10  km.  The  apparent  velocities  of  the  computed  wave  indicate  that  they  begin  seeing  the 
faster,  deeper  layer  at  distances  between  about  50  and  70  km. 

Noise  was  not  added  to  the  seismograms,  but  it  is  clear  that  P-waves  in  the  low-Q 
model  would  be  difficult  to  detect  in  the  presence  of  noise  at  all  three  frequencies  at 
distances  greater  than  about  150  kin  for  small  events  and  that  the  onsets  are  so  emergent 
that  they  would  not  be  useful  for  precise  event  location  even  at  much  shorter  distances. 
The  high-Q  model  produces  records  in  which  amplitudes  attenuate  much  more  slowly 


and  could,  at  the  higher  frequencies,  probably  be  detected  at  distances  beyond  200  km  if 
the  noise  level  is  not  high.  The  sharpness  of  the  wave  onsets  for  the  high-Q  model  also 
makes  those  waves  more  useful  for  event  location. 

Figure  4  plots  vertical-component  P-wave  amplitudes  as  a  function  of  distance  for 
the  three  frequencies  for  both  our  high-  and  low-Q  models.  For  each  frequency  the 
amplitude  maxima  for  the  low-  and  high-Q  models  occurs  at  about  the  same  distance, 
where  that  distance  increases  with  decreasing  frequency.  Amplitude  variations  with 
distance  show  that  1-Hz  waves  have  the  largest  amplitudes  of  the  three  frequencies  at 
distances  begin  at  about  70  km  for  the  high-Q  model  and  at  about  65  km  for  the  low-Q 
model.  Amplitude  for  the  1-Hz  waves  are  twice  or  more  as  large  as  for  the  other 
frequencies  at  all  greater  distances.  At  distances  of  40-50  km  the  amplitudes  of  the  5-Hz 
waves  are  more  than  four  times  greater  than  the  1-Hz  waves  and  about  1.5  times  greater 
than  the  3-Hz  waves. 

Similar  distributions  were  obtained  for  the  radial  component  of  computed  P-waves 
(Figure  5).  Peak  amplitudes  of  radial-component  waves  are  higher  than  for  the  vertical- 
component,  but  this  is  only  because  of  the  constant  velocity  in  the  upper  layer  of  the 
model.  With  a  more  realistic  model  in  which  velocities  increase  with  depth  the  vertical- 
component  amplitudes  would  become  much  larger  than  those  obtained  in  this  preliminary 
work. 


Conclusions 

Computations  of  vertical-  and  horizontal-component  ground  motion  from  a  source  at 
5  km  depth  and  an  assumed  value  of  0.5  for  the  frequency  dependency  parameter  in  the 
upper  crust  have  allowed  us  to  determine  the  variation  of  P-wave  amplitudes  with 
distance  in  a  simple  two-layer  model.  Synthetic  seismograms  indicate  that  it  may  be 
difficult  to  detect  small  events  at  larger  distances,  especially  if  noise  is  present,  for  the 
low-Q  model  and  that  emergent  waveforms,  especially  at  1  Hz,  will  degrade  attempts  to 
locate  the  events  over  a  large  distance  range.  At  higher  frequencies  (3  and  5  Hz)  the 
onsets  are  sharper,  thus  yielding  better  locations,  but  the  waves  attenuate  very  rapidly  for 
the  low-Q  relative  to  the  high-Q  model.  Synthetics  at  the  higher  frequencies  decay  very 
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rapidly  for  the  low-Q  model  relative  to  the  high-Q  model  whereas  synthetics  at  1  Hz 
show  less  difference  for  the  model  we  have  tested. 

Amplitude  maxima  shift  to  larger  distances  with  decrease  in  frequency.  For  our 
selected  model  amplitude  maxima  occur  at  distances  of  45,  60  and  75  km  for  frequencies 
of  5,  3  and  1  Hz,  respectively.  Plotted  amplitudes  for  the  higher  frequencies  are  large  at 
distances  near  and  inside  the  critical  distance  but  fall  off  rapidly  at  distances  greater  than 
about  70  km.  Selection  of  other  values  for  the  frequency  dependence  parameter  are  likely 
to  produce  far  difference  results  for  amplitude  attenuation  with  distance  for  the  high- 
frequency  waves  of  this  study.  With  our  selected  Rvalue  maximum  amplitudes  at  5  Hz 
are  4  times  and  1.5  times  greater  than  at  1  Hz  and  3  Hz,  respectively. 
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Table  Caption 

Table  1.  Velocity-Q  models  utilized  for  the  computations  of  this  study. 


Figure  Captions 

Figure  1.  Vertical-component  seismograms  bandpass  filtered  around  1  Hz. 

Figure  2.  Vertical-component  seismograms  bandpass  filtered  around  3  FIz. 

Figure  3.  Vertical-component  seismograms  bandpass  filtered  around  5  Hz. 

Figure  4.  Amplitude  variation  with  distance  of  P-wave  vertical-component  ground 
motion  at  frequencies  of  1,  3  and  5  FIz  for  high-  and  low-Q  models. 

Figure  5.  Amplitude  variation  with  distance  of  P-wave  radial-component  ground 
motion  at  frequencies  of  1,  3  and  5  Hz  for  high-  and  low-Q  models. 
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Table  1 


Layer  1 

Layer  2 

Thickness  (km) 

10 

CO 

P-wave  vel.  (km/s) 

6.1 

6.7 

S-wave  vel.  (km/s) 

3.1 

3.7 

Density  (gm/cc) 

3.5 

3.8 

QP 

Low-Q  model 

100 

High-Q  model 

1000 

600 

Qs 

Low-Q  model 

50 

70 

High-Q  model 

500 

300 

Epicentral  Distance  [km] 


Epicentral  Distance  [km] 


Epicentral  Distance  [km] 


Amplitud  [*10'9] 


Amplitud  [*10‘8] 
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